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ty 

David Fisher 
and 

Theodore L. Felsentrager 


SUMMARY 



It is shown that the long range variations of the orbital ele- 
ments of a lunar satellite in libration is periodic, but non- analytic, 
in nature. A method for finding the period of the motion, as well 
as the variations of the elements with time, is described. 
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THE LIBRATION OF A LUNAR SATELLITE 


INTRODUCTION 

The basic equations for the long range motion of a lunar satellite lead to 
two types of motion: circulation and libration. These equations involve the 
argument of pericenter and the eccentricity (or some function of the eccentricity) 
as the variables (see Reference 1). When circulation occurs, all values of the 
argument of pericenter ranging between 0 and 277, as well as multiples of this 
range, are permitted. However, when the range of values of g is restricted to 
a limited set, the motion is called libration. These two types of motion are il- 
lustrated by a pendulum which is free to rotate about a horizontal axis. When the 
pendulum forms complete revolutions about this axis, the motion is circulation; 
however, if the pendulum oscillates about a vertical axis (e.g., the pendulum of 
a clock), the motion is called libration. 

The equations describing the long period motion of a lunar satellite are 
similar to the equations for the stellar problem of three bodies discussed by 
E. W. Brown (Reference 2) over thirty years ago. In an interesting report, 
Frances Frost (Reference 3) has shown that the method of harmonic analysis 
used by Brown is well adapted to the calculation of the motion of lunar satellites 
in the case of circulation. However, this method cannot be applied when libration 
takes place. 

It is possible to determine the argument of pericenter g as a function of 
its time derivative g from the energy equation and the equation for g given 
below. It has been shown in treatises on non-linear mechanics (e.g., see Ref- 
erence 4) that in the (g,g) or displacement-velocity space, libration leads to a 
closed curve or cycle, while circulation leads to an open curve. 

It will be shown below that the cycle described by g consists of two distinct 
branches. This introduces an unexpected asymmetry in the results, leading to a 
discontinuity in the derivatives of the cycle where the two branches join. 

The argument of pericenter g , and consequently g and e , can be found as 
functions of time by integration over these branches. The period of oscillation 
of the other elements of the lunar satellite may be found by integrating over the 
entire region of the cycle in the displacement-velocity space. 
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In order to make the analysis concrete, a specific example is chosen in 
this report and the calculations necessary to obtain the variations of all the ele- 
ments as functions of time are given in detail. 


ANALYSIS IN THE PHASE SPACE OF A PENDULUM 

In order to provide motivation for the analysis of a lunar satellite in the 
displacement-velocity space (or phase space, as it is often called), we first 
examine the linear motion of a pendulum. We assume that the angle of oscilla- 
tion of the pendulum about the vertical is small enough so that the equation of 
motion can be written in the form 


x + a 2 x = 0, 


( 1 ) 


where x is the angle of displacement from the vertical and x is the second 
derivative with respect to time. 

The solution to Equation (1) can be written in the form 

x = ± b cos (a t + c) 

x = + ab sin (at + c) , 


( 2 ) 


where b and c are constants of integration and a is the frequency of oscillation. 
Equations (2) also satisfy the relation. 


x 2 


b 2 a 2 b 2 


= 1 . 


( 3 ) 


Equation (3) is an integral of the motion, and is often called the energy integral. 

The analysis in the displacement -velocity space begins with Equation (3) . 

We illustrate the method followed in the more complicated lunar satellite problem 
by using Equation (3) to derive the results given by Equations (2). 

Equation (3) represents an ellipse in the x , x space (see Figure 1); the 
semi-major and semi-minor axes are b and ab, respectively. The motion 
around the ellipse is taken positive in a counterclockwise direction, and the con- 
vention for the four quadrants is taken as shown in Figure 1. 
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In addition, the indefinite integral from Equation ( 5 ) gives 

x = ab sin (at + c) , 

where c is a constant of integration. Finally, 

x = - b cos (at + c) , 

where the negative sign is the result of the sign convention chosen in Equation (4). 

This example illustrates the basic principle adopted in this report for the 
analysis of the libration of a lunar satellite. 


THE ENERGY INTEGRAL AND THE EQUATION OF 
MOTION OF THE ARGUMENT OF PERICENTER 

The expressions for the energy integral and the equation of motion of the 
argument of pericenter may be found in References 5 and 6. The energy equa- 
tion can be written in the form 

f i i 7 !) + V 3 f 2 (V) sin 2 g = 0; (7) 


the variables are 77 and g , where 77 = /l - e 2 ( e is the eccentricity) , and g is 
the argument of pericenter. The functions f x (77) and f 2 (77) are defined by 

f, (77) = - 12— t] 7 + A (10 + t 

1 q q 

f (77) = 30 — (t] 2 - 1) (77 2 - v 2 ) 

2 q 

In addition, the equation for the motion of the argument of pericenter is 



The quantities appearing in Equations (8) and ( 9 ) are defined in Appendix A. 


. 2n 4 C 

5 V 1 ) 

3 nL 


77 s - 2 a 2 t } 2 + 6 a 2 v 2 


( 8 ) 
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THE CYCLE IN PHASE SPACE 


We have indicated above that libration occurs when the argument of peri- 
center g is restricted to a finite range of values such that a plot of g versus g 
in the phase space (called a cycle) is a closed curve (e.g., Figure 1). 

In order to illustrate such a case of libration, we shall examine Equations 
(7), (8), and (9) for the initial conditions 

e = .1, cosi = .71066905, g = 90°, a = 7 4822577 moon radii. 


For these initial conditions, Equations (8) become 

fj (V) = -22.947701 t? 7 + 22-551912 v 5 - -8609575 v 2 + 1-29143625 

( 10 ) 

f 2 (rj) = 22.684626 (2v 2 - 1) (V 2 - !)• 


For the convenience of the reader, we list the values for the following constants: 
— = 1.9123084 x 10" s , a, = .43047875 x 10‘ 5 

q 





n = 3.8097588 x 10“ 2 conical (Vanguard) units. 


Substitution of Equations (10) into Equation (7) provides a relation between 
the eccentricity e and the argument of pericenter g . This functional relation is 
plotted in Figure 2, which indicates a distorted "ellipse." The variable g 
ranges between the values 75?96481 and 104?03519, while e varies between .1 
and .30694755. The major axis of this "ellipse" is at e = .17962767, which may 
be considered a "mean" value of sorts. Adopting the quadrant convention shown 
in Figure 1, the numerical values of e and 90° — g are listed in Table 1. 

The next step is to find g as a function of g , which can be readily accom- 
plished by substituting the values of e and 90° -g in Table 1 into Equation (9). 

It is found that the values of g in 2.8573191 x 10 ” 7 canonical units vary from 
—1.0396376 to 1.0950486— these values are also listed in Table 1. . 

Figure 3 shows the 90°— g, g cycle derived from the values in Table 1. It 
is seen that the figure departs only slightly from an ellipse, with the branch 
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90°-g (deg) 


Figure 2-e vs. 90° — g. 

consisting of quadrants I and II differing a little from the branch contained in 
quadrants III and IV. 


EXPANSIONS IN PHASE SPACE 


The results illustrated by Figure 3 may be utilized to find analytic expres- 
sions for the cycle, and subsequently to express the elements of the lunar satel- 
lite as functions of time. 


It is convenient to define the non-dimensional variables u , , and w 2 by 


- 90° -g w g 

14*03519 1 2.9705764 x 10' 7 


w„ 


g 

3.1289033 x 10" 7 


6 






90° -g (deg) 

Figure 3— — -gx 10 5 vs. 90° - g 
3n 

where g is in degrees and g is in canonical units. Values for u, Wj, and w 
appear in Table 1. The angle 14?03519 is the variation of g about 90°, while 
2.9705764 x 10" 7 and —3.1289033 x 10" 7 are, respectively, the maximum and 
minimum values of g . Using the fact that the plot in Figure 3 consists of two 
branches, the following expansions were derived by least squares from the 
values in Table 1: 


u = - {/l - w 7 - .22131760 w 7 + .20999084 w« } for quadrant I 
u = ^ /l - w 7 - -22131760w 7 + .20999084 w* ^ for quadrant II 
u = {/ 1 - w 7 - .26127200 w 7 + .24764299 jfor quadrant III 
u = - { - w 7 - -26127200 w 7 + .24764299 w* } for quadrant IV (H) 
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The expansions in Equations (11) were obtained by means of a desk calculator, 
and are therefore carried out to only a limited number of terms. However, they 
do represent the cycle within an accuracy of about 1%. 

The next step in the process is to utilize Equations (11) to find the periodicities 
of the system, as well as expansions in terms of time. 


EXPANSIONS OF g, g, AND e IN TERMS OF TIME 
The differential equation for the time is 

( 12 ) 

g 

integration is performed over the cycle of Figure 3, with the aid of Equations (11). 

Since there are essentially two different expansions in Equations (11) (one for 
each branch), we can determine two different periodicities in g and g. Thus, 
we define 



where the integration is performed over quadrant I or II, and 



(13) 


(14) 


where the integration is performed over quadrant III or IV. Then, by substituting 
Equations (11) into Equations (13) and (14), we find that 

Tj = .82462207 x 10 6 [sin _1 Wj + .44263520 - .27998779 w3] 

= 1.4294360 x 10 6 canonical time units / 
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T 2 = .78289502 x 10 6 [sin -1 w 2 + .55254400 w 2 - .33019065 ] 1 

= 1.3803611 x 10 6 canonical time units. 

It follows that the period taken over the entire cycle of Figure 3 is 

T = 2 (Tj +T 2 ). 

Equations (11) and (12) can be utilized to obtain the time histories of g and 
g. Thus, we find that 

Ij = x .82462207 x 10 6 [ sin' 1 w x + .44263520w x - ,27998779w f ] +const. 

T i 

77 (15) 

= •=“ (t - 1 0 ) + const. 

A i 


represents w x as a function of time in quadrant I, and - g x represents w l in 
quadrant II. Similarly, 

g 2 = — x .78289502 x 10 6 [ sin' 1 w 2 + .52254400w 2 - .33019065wS j +const. 

t 2 

(16) 

= ~ (t - t Q ) + const. 

l 2 


represents w 2 in quadrant m, and - g 2 represents w 2 in quadrant IV. 

Now, let us choose t - t Q to be zero at the time of the initial conditions (i.e., 
the beginning of quadrant II); in addition, let g x and g 2 increase with time with 
mean motions 77 /T x and 77 /T 2 , respectively. Further, let g x = 90° at t - t Q = 0 . 
Values of g x , g 2 , and t - t Q appear in Table II. 

Figures 4 and 5 are plots of u and e, respectively, vs. t - t Q . Least squares 
analyses yield the following results: 


Uj ~ — 1.02099718 cos g x + .01610990 cos 3 g x + -00649749 cos SgjvlT) 
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80 



0 32 



li 


Figure 5— e vs. 


e = .17962767 - .07962767 (1.03054138 sin g t + .03349887 sin 3 g 1 

+ .00518834 sin 5 g 2 ) (18) 

for quadrants I and II, and 

u 2 = - 1.02331813 cos g 2 + .01975752 cos 3 g 2 + .00608506 cos 5 g 2 (19) 

e = .17962767 + .12731988 (- .95873486 sin g 2 + .04685613 sin 3 g 2 
+ .00286600 sin 5 g 2 ) (20) 

for quadrants III and IV. Equations (17), (18), (19), and (20) represent the data 
to within about 2% . 


I AND h 

Equations for the motions of the mean anomaly Z and the Delaunay variable 
h may be obtained from References 5 and 6 — they are as follows: 



10 + 3*’ - 6v 2 - 15 (2-V)(V->' : ) , in .; 


V‘ 


h = 




3 

4 q 


0 10 (1 - T7 2 ) . 2 

2 + — — sm^ g 



la 2 77 - 3 ^-\ 

4 2 J 

( 21 ) 

( 22 ) 


(See Appendix A for definitions of constants and variables.) 

Values for A \ /n x io 5 and Ah/nvx 10 s appear in Table in ( = ^- n , 
and Ah = h + n^). i n addition, Figures 6 and 7 show plots of A^/ n x 10 5 and 
Ah/nv x 10 s , respectively, vs. time. These can be integrated numerically to 
find the perturbations in i and h . 
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3.7 





CONCLUSIONS 


The type of motion known as libration is possible from lunar satellites — by 
examining the energy equation and the equation of motion of the argument of 
pericenter, initial conditions leading to libration can be found. 

The methods of non-linear mechanics are useful for investigating the details 
of the motion. The procedure chosen for the example in this report may be 
adapted to electronic machine calculations. 

The results presented here differ from those obtained through the use of 
usual perturbation theories in the following ways: 

1 . The variation in e is much larger than the order of magnitude of the 
disturbing function, which is about 10 " s . 

2. The variation in the argument of pericenter can be explained without a 
secular term. 

Finally, the fact that substantial changes may occur in the eccentricity of a 
lunar satellite in libratory motion should be taken into account in the planning of 
various types of orbits. 
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Appendix A 
SYMBOLS 

a = semi-major axis of satellite's orbit 
e = eccentricity of satellite's orbit 

i = inclination of satellite's orbital plane to moon's equatorial plane 
■i = mean anomaly of satellite 
g = argument of pericenter of satellite's orbit 
h = 

Q = longitude of ascending node of satellite's orbit 
v e = mean longitude of earth, measured on moon's equator 
n = mean motion of satellite 

(J. = n 2 a 3 
L = /JTa 
v = /l - e 2 cos i 
t] = /l - e 2 

a fi = semi -major axis of moon's orbit - 384,400 Km. 

= mean radius of moon - 1738 Km. 

J 2 = second zonal harmonic coefficient of moon - 2.41 x io - 4 
n! = mean motion of X 

© 

- mass of moon 

q = 1 + 

mass of earth 
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Table I 


Quadrants I and II 


90° -g (deg.) 

e 

-f i x i° s 

3n 

u 

0 

.10000000 

-1.0396376 

0 

T 1.40225 

.10027263 

-1.0328040 

t.09990958 

T 2.48788 

.10086620 

-1.0180892 

t.17726016 

T 5.08020 

.10378266 

- .94899765 

t.36196161 

T 7.66137 

.10938007 

- .82966344 

T. 54586864 

T 9.98252 

.11811435 

- .67139774 

T. 71124937 

Til. 83 747 

.13010669 

- .49435338 

T. 84341 359 

T12. 85549 

.14106736 

- .36056622 

T. 91 594699 

T13. 79363 

.16072025 

- .16220364 

t.98278898 

Tl4. 03390 

.17816846 

- .01190701 

t.99990809 

t14. 03519 

.17962767 

0 

Tl 


Quadrants m and IV 


90° -g (deg.) 

e 

-ig x 105 
3n 

u 

±14.03519 

.17962767 

0 

±1 

±14.02577 

.18360556 

.03209762 

±.99932883 

±13.57598 

.20860489 

.22663012 

±.96728153 

±12.65501 

.23080511 

.39860814 

±.90166289 

±11.92278 

.24310492 

.49708257 

±.84949189 

±10.39403 

.26224988 

.65845333 

±.74056924 

± 8.36212 

.28000000 

.82007223 

±.59579671 

± 7.33139 

.28676820 

.88544486 

±.52235733 

± 4.47007 

.29980660 

1.0181686 

±.31849017 

± 2.61962 

.30453859 

1.0687611 

±.18664656 

± 1.56060 

.30609802 

1.0857364 

±.11119194 

0 

.30694755 

1.0950486 

0 


-1 

-.99342694 

-.97927316 

-.91281582 

-.79803139 

-.64579979 

-.47550549 

-.34681914 

-.15601941 

-.01145304 

0 


0 

.02931159 

.20695896 

.36400954 

.45393654 

.60130056 

.74889117 

.80858956 

.92979307 

.97599422 

.99149609 

1 
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Table II 


t-t 0 (10 4 days) 


Quadrant II 

Quadrant I 

0 

5.2477884 

.0863580 

5.1614304 

.1509922 

5.0967962 

.3020400 

4.9457484 

.4608704 

4.7869180 

.6321681 

4.6156203 

.8143666 

4.4334218 

.9529006 

4.3218878 

1.1618342 

4.0859342 

1.3221378 

3.9256506 

1.3348611 

3.9129273 


t-t 0 (io 4 

days) 

Quadrant III 

Quadrant IV 

1.3348611 

3.9129273 

1.3674856 

3.8803028 

1.5641943 

3.6835941 

1.7346646 

3.5131238 

1.8301755 

3.4176129 

1.9837445 

3.2640439 

2.1383611 

3.1094273 

2.2046564 

3.0431320 

2.3688471 

2.8789413 

2.4711777 

2.7766107 

2.5313346 

2.7164538 

2.6238942 

2.6238942 


g x (deg.) 


Quadrant n 

Quadrant I 

90 

90 

95.822496 

84.177504 

100.180309 

79.819691 

110.364367 

69.635633 

121.073147 

58.926853 

132.622511 

47.377489 

144.906837 

35.093163 

154.247180 

25.752820 

168.334056 

11.665944 

179.142166 

0.857834 

180 

0 

g 2 

(deg.) 

Quadrant III 

Quadrant IV 

180 

0 

182.277837 

357.722163 

196.011995 

343.988005 

207.914188 

333.085812 

214.582740 

325.417260 

225.304893 

314.695107 

236.100188 

303.899812 

240.728914 

299.271086 

252.192670 

287.807330 

259.337369 

280.662631 

263.537508 

276.462492 

270 

270 
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Table in 


Quadrants I and n 


A£ x 10 s 

10 s 

n 

n v 

2.0223381 

-3.6754823 

2.0177568 

-3.6762854 

2.0079104 

-3.6780420 

1.9620501 

-3.6868374 

1.8845147 

-3.7044884 

1.7859026 

-3.7340862 

1.6835939 

-3.7789088 

1.6141463 

-3.8242497 

1.5284436 

-3.9165212 

1.4805491 

-4.0109717 

1.4774378 

-4.0194326 

aadrants in and IV 

1.4774378 

-4.0194366 

1.4695560 

-4.0429523 

1.4368206 

-4.2067104 

1.4264079 

-4.3771699 

1.4264217 

-4.4828268 

1.4332899 

-4.6648510 

1.4462603 

-4.8546807 

1.4527512 

-4.9329060 

1.4675579 

-5.0934689 

1.4736671 

-5.1551381 

1.4757658 

-5.1758757 

1.4769324 

-5.1872620 
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